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We present an update of our systematic analysis of all Gamma-Ray Burst (GRB) afterglow data, now published 
through the end of 2004, in an attempt to detect the predicted supernova light component. We fit the observed 
photometric light curves as the sum of an afterglow, an underlying host galaxy, and a supernova component. 
The latter is modeled using published UBVRI light curves of SN 1998bw as a template. The total sample 
of afterglows with established redshifts contains now 29 bursts (GRB 970228 - GRB 041006). For 13 of them 
a weak supernova excess (scaled to SN 1998bw) was found. In agreement with our earlier result E2 we find 
that also in the updated sample all bursts with redshift <0.7 show a supernova excess in their afterglow light 
curves. The general lack of a detection of a supernova component at larger redshifts can be explained with 
selection effec ts. These results strongly support our previous conclusion based on all afterglow data of the years 
1997 to 2002 |47l that in fact all afterglows of long-duration GRBs contain light from an associated supernova. 


1. Introduction 

Significant progress towards understanding the na¬ 
ture of GRBs and their progenitors came with the 
discovery of GRB afterglows in 1997 IS 23 . First 
observational evidence for the underlying source pop¬ 
ulation was provided by GRB 970828, which showed 
a bright X-ray afterglow but no optical counterpart 
down to faint magnitudes 0 - This led to the sug¬ 
gestion that the optical light was blocked by cosmic 
dust in the GRB host galaxy, linking the burster to 
a dusty star-forming region, i.e., most likely to the 
explosion of a massive star [0. The discovery of a 
near-by type Ibc supernova (SN 1998bw) in the error 
circle of the X-ray afterglow for GRB 980425 0,01 > 
provided strong support for this idea, and is consis¬ 
tent with our current understanding of type Ibc SNe 
and their progenitors (e.g., 110 ) 

From the observational site, the supernova picture 
is further supported by the fact that all GRB hosts 
are star-forming, and in some cases even star-bursting 
galaxies (e.g., 0, 4(|). Evidence for host extinction by 
cosmic dust in GRB afterglows and the discovery of 
an ensemble of optically ’dark bursts’ (for a recent dis¬ 
cussion, see if, ( 23 , 0 ) also is consistent with the pic¬ 
ture that GRB progenitors are young, massive stars. 
Furthermore, for several GRB afterglows X-ray lines 
may hint at a period of nucleosynthesis preceding or 
accompanying the burst yj, 0, (3C|. The positions of 
the afterglows with respect to their hosts also favors 
a relation to young, massive stars to GRBs [|j ■ 

As a natural consequence of a physical relation be¬ 
tween the explosion of massive stars and GRBs super¬ 
nova light should contribute to the afterglow flux, and 
even dominate under favorable conditions. The most 
convincing example is GRB 030329 0 at 2=0.1685 
0 with spectral confirmation of supernova light in its 
afterglow 0, 123, 0, 0] ■ Spectroscopic evidence for 
SN light in a GRB afterglow was later also reported 
for GRB 021211 @], GRB 031203 0 and most re¬ 



Figure 1: The late-time bump in the optical afterglow of 
GRB 970228, the first optical afterglow ever found (data 
collected from the literature). Interpreted as a signature 
from a SN explosion makes this event one of the most 
distant core-collapse SNe ever seen at that time. See also 

0IH. 

cently for XRF 020903 0. 

In contrast to direct spectroscopic evidence, several 
cases of photometric indication of extra light in GRB 
afterglows have been reported, starting with the pi¬ 
oneering work on GRB 980326 JJ. Inspired by this 
finding, the discovery of extra light in archived data 
of the afterglow of GRB 970228 0,|23 : ] made it clear 
that a search for late-time bumps in optical after¬ 
glow light curves provides a powerful tool to constrain 
or even reveal the nature of the underlying sources. 
Since then various groups successfully fit SN 1998bw 
templates to explain these late-time bumps (e.g., i), 
the most convincing case being that of GRB 011121 

i00- 

The goal of our study is to search for supernova 
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bumps in GRB afterglow light curves using a sys¬ 
tematic approach, allowing us to draw statistically 
founded conclusions on the physical properties of this 
new class of GRB-SNe in particular and on the GRB 
progenitors in general. We collected from the liter¬ 
ature all available photometric data on GRB after¬ 
glows (including our own data), checked them for pho¬ 
tometric consistency, and re-analyzed the data in a 
consistent manner. Here we report on the status of 
our study for all bursts that occurred by the end of 
2003, supplementing and expanding our previous re¬ 
sults (H 3 , in the following paper I). 



2. Numerical approach 

We model the light curve of the optical transient 
(OT) following a GRB as a composite of afterglow 
(AG) light, supernova (SN) light, and constant light 
from the underlying host galaxy. The flux density, F u , 
at a frequency v is then given by 

F GT (t) = F^ G (t) + k F„ SN (t/s) + E„ host . (1) 

Here, the parameter k describes the observed bright¬ 
ness ratio (in the host frame, i.e., including the cos¬ 
mological /v-correction) between the GRB-supernova, 
and the SN template (SN 1998bw) in the considered 
photometric band (in the observer frame). We allowed 
k to be different in every photometric band, but within 
a band independent of frequency. The parameter s is 
a stretch factor with respect to the used template. We 
have also explored the consequences of a shift in time 
between the onset of the burst and the onset of the 
supernova explosion, as implied by some theoretical 
models [Sj. Then, in Eq. m F™(t/ S ) was replaced 
by F® N (t + r). Here, t = 0 refers to GRB 980425/SN 
1998bw (2l| . If r < 0 the SN preceded the onset of 
the GRB. 

Following p) and jftfj , we describe the afterglow 
light curve by a broken power-law, 

F^ G (t) = const [{t/t b ) ai n + (t/t b ) a * n }- 1 ' n , (2) 

with const=2 1 / n F^ G {t b ). Here t is the time after the 
burst (in the observer frame), aq is the pre-break de¬ 
cay slope of the afterglow light curve, a .2 is the post¬ 
break decay slope, and t b is the break time. The pa¬ 
rameter n characterizes the sharpness of the break; a 
larger n implies a sharper break. If no break is seen 
in the data then aq = ol-i and Eq. © simplifies corre¬ 
spondingly. 

The results of this numerical procedure were com¬ 
pared with corresponding results published by oth¬ 
ers 13 S, and we found close agreement. We used 
this procedure to predict the color evolution of GRB 
030329/SN 2002dlr |4(|, and obtained a very good nu¬ 
merical fit for the light curves of GRB-SN 011121 (l5j . 


Figure 2: The afterglow of GRB 011121 showed a very 
clear signature of a late-time bump rising some days after 
the burst. Shown here are data obtained with the 
telescopes at ESO, Chile, and with the Hubble Space 
Telescope Q). The bump can be modeled well by an 
underlying SN component at the redshift of the burster 
(^=0.36) with a peak luminosity of about 80% of the 
peak luminosity of SN 1998bw. Note that in the figure 
the flux from the underlying host galaxy was subtracted 
from the data (for details, see El). 



t (days after the burst) 
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10 2 


io' 





Figure 3: Sketch of the hidden SN bump in the afterglow 
of GRB 030329. Various re-brightening episodes of the 
genuine afterglow, in combination with a relatively late 
break-time of the light curve, made the photometric 
signature for the underlying SN explosion very small. 
Presumably, with no spectroscopic evidence at hand, the 
SN had easily been missed in the data. For the light 
curve fit 0:2 was fixed at 2.5. 


The limitations of the procedure are given by the cho¬ 
sen photometric band in combination with the redshift 
of the burster. Once we can no longer interpolate in 
between the UBVRI bands, but have to extrapolate 
into the UV domain (cf. j|), results become less ac¬ 
curate. For more details see paper I. 

Before performing a numerical fit, the observational 
data was corrected for Galactic extinction along the 
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Table I The input sample of GRB afterglows. Redshifts 
were taken from the literature. 


GRB 2 

GRB 2 

GRB 2 

GRB 2 

970228 0.695 

970508 0.835 

971214 3.42 

980703 0.966 

990123 1.600 

990510 1.619 

990712 0.434 

991208 0.706 

991216 1.02 

000301C 2.04 

000418 1.118 

000911 1.058 

000926 2.066 

010222 1.477 

010921 0.450 

011121 0.362 

011211 2.140 

020405 0.69 

020813 1.25 

020903 0.251 

021004 2.3 

021211 1.01 

030226 1.986 

030323 3.372 

030328 1.520 

030329 0.169 

030429 2.658 

031203 0.106 

041006 0.716 


line of sight using the COBE maps _37]. This also 
holds for SN 1998bw, where we assumed E(B — V)= 
0.06 mag. We calculated the Galactic visual extinc¬ 
tion according to Ay al = 3.1 E(B — V ), whereas the 
extinction in U and B were obtained via [36j], and in 
R c and I c by means of the numerical functions com¬ 
piled in 

Most of the light curves we investigated have been 
followed in more than one photometric band. For 
each of these GRBs we chose the best-sampled light 
curve as a reference light curve for the fit in the 
other photometric bands. In all cases this was the 
R band light curve. We always assumed that after¬ 
glows are achromatic, in reasonable agreement with 
observational data (e.g., |30|)- For every individ¬ 
ual GRB, the afterglow parameters ai,a2,ft, and n 
(Eq.H are then the same for all photometric bands. 
Consequently, once we fit the reference light curve of 
an optical transient and deduced the corresponding 
afterglow parameters, we treated them as fixed pa¬ 
rameters when fitting the light curves of the optical 
transient in other photometric bands. In the fit the 
degrees of freedom are reduced correspondingly. 


3. Results and Discussion 

The input sample consists of 29 bursts (GRB 
970228 - GRB 041006) with established redshifts and 
good enough photometric data in order to search for 
a late-tinre bump in their afterglows (Table 0 > with 
the most recent data for GRB 041006 |42ll . These are 
eight bursts more than in our previous study for all 
bursts observed by the end of 2002 (paper I). Among 
these 29 bursts are 13 for which a late-time bump 
was found. This includes now also XRF 020903 (as 
already noted in l38ll and now spectroscopically con¬ 
firmed J3fj)- Note that the requirement of a known 
redshift excludes GRB 980326 as well as XRF 030723 
from this list, which both showed a strong late-time 
bump. On the other hand, as already noted in Jj, 


one can in principle constrain the redshift of a burster 
by fitting a redshifted SN component to the observed 
late-time bump in its afterglow light curve. 
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Figure 4: 

The deduced peak luminosities of all GRB-SNe 


in units of the peak luminosity of SN 1998bw. All data 
refer to the R c band (in the observer frame). The dotted 
line corresponds to SN1998bw. The parameter Am 
equals -2.5 log k, which measures the magnitude 
difference at maximum light between the GRB-SN and 
SN 1998bw in the corresponding wavelength regime. 

Note that the data are not corrected for a possible 
extinction in the GRB host galaxies. GRB 980703 is not 
included here and in the following figures because in this 
case the stretch factor was not allowed to vary freely. 

Again, our key finding is photometric evidence of 
a late-time bump in all GRB afterglows with a red¬ 
shift 2 <, 0.7. We interpret this bump as light from 
an underlying supernova, and model this component 
as a redshifted version of SN 1998bw. The deduced 
luminosities for these GRB-SNe (not including ex¬ 
tinction corrections for the host galaxy), normalized 
to SN 1998bw are listed in Table [H] The width of 
the distribution of the SN peak luminosities (in units 
of the peak luminosity of SN 1998bw) spans over 2 
photometric magnitudes with a pronounced maximum 
around k = 0.5...1 (Figs. 0J0. The potential SN re¬ 
lated to the X-ray flash 020903 is not unusual with re¬ 
spect to its peak luminosity. Interestingly, SN 1998bw 
is at the bright end of the GRB-SNe distribution (as 
already noted in paper I). Only the SNe related to 
GRBs 030329 and 031203 might have been slightly 
more luminous at peak brightness. No correlation 
was found of the deduced SN luminosities with the 
redshift or any afterglow parameter. Note, however, 
that we cannot exclude the existence of such a corre¬ 
lation since in most cases when a SN was found there 
is a lack of early time data in the optical light curve 
(resulting in an unknown break time ft and, hence, an 
unknown parameter oq; Eq. n. 

Figure IH1 shows the distribution of the correspond- 
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Table II Best-fit parameters for the SN component found in GRB afterglows with known redshift. Columns: (1) and 
(2) GRB and redshift; (3) photometric band, in which the light curve was fitted; (4) central wavelength of the 
photometric band in the host frame in units of nm, adopting a wavelength of 659 nm for the R c band; (5) peak 
luminosity of the fitted SN component in the corresponding wavelength band (observer frame) in units of SN 1998bw, 
after correction for Galactic extinction; (6) stretch factor s (Eq.Q; (7) goodness of fit per degree of freedom; (8) and 
(9) the same as (5) and (7) for s = 1. The low y 2 /d.o.f. for GRB 970228 is due to the small number of data points. 
Note that we reduced all data with our own numerical procedure, so that slight differences to the results obtained by 
others do naturally exist. 


GRB 

z 

band 

^host 

k 

s 

Xd.o.f. 

k if ,s=l 

Xd.o.f. 

970228 

0.695 

Rc 

389 

0.40T0.24 

1.46±0.80 

0.70 

0.33T0.30 

0.71 

980703 

0.966 

Rc 

335 

- 

- 

- 

1.66±1.22 

0.78 

990712 

0.434 

Rc 

459 

0.48±0.10 

0.89±0.10 

1.00 

0.43±0.08 

1.01 

991208 

0.706 

Rc 

386 

0.90±0.35 

1.12±0.28 

1.64 

1.02±0.32 

1.56 

000911 

1.058 

Rc 

320 

0.87T0.39 

1.49±0.33 

0.75 

0.51T0.43 

1.14 

010921 

0.450 

Rc 

454 

0.68T0.48 

0.68T0.28 

0.42 

0.43±0.10 

0.78 

011121 

0.360 

Rc 

484 

0.79T0.06 

0.85±0.06 

0.92 

0.74T0.05 

1.32 

020405 

0.695 

Rc 

389 

0.74±0.17 

0.98T0.17 

5.26 

0.72±0.11 

4.86 

020903 

0.251 

Rc 

527 

0.46±0.44 

1.47±0.88 

1.52 

0.37T0.41 

1.22 

021211 

1.006 

Rc 

328 

0.97±0.87 

0.74±0.23 

2.68 

0.52±0.34 

2.65 

030329 

0.169 

Rc 

563 

1.13±0.33 

0.82±0.13 

3.10 

0.98±0.01 

4.49 

031203 

0.106 

Rc 

596 

1.65±0.41 

1.14±0.16 

0.04 

1.75±0.19 

0.24 

041006 

0.716 

Rc 

384 

0.91±0.05 

1.38±0.06 

1.27 

1.21±0.07 

1.95 



o.o 

o. 


991208 041006 


0.5 1.0 1.5 

s (stretch factor ) 


Figure 5: The distribution of the peak brightness of all 
GRB-SNe in units of the corresponding peak brightness 
of SN 1998bw in the R c band (in the observer frame). 
The corresponding stretch factor is shown in Fig.|S| 


ing stretch factor s. Since no fit was possible for the 
afterglow of GRB 980703 with s being a free parame¬ 
ter, this burst is not included in Figs. l4IOl The mean 
value of s is 1.0, i.e., identical to SN 1998bw. 

Instead of introducing a stretch factor to have more 
freedom in the variety of GRB-SNe, one can also fol- 


Figure 6: The distribution of the corresponding 
parameter s (Eq. 0 describing a stretching of the SN 
light curve relative to those of SN 1998bw (for which by 
definition s = 1, dotted line). In general, s < 0 (s > 0) 
means that the evolution of the light curve of the SN was 
slower (faster) than those of the light curve of SN 1998bw 
in the corresponding wavelength band. The mean value 
is s = 1.0. 

low 44] and search for evidence of a time delay be¬ 
tween the burst and the SN. According to this model, 
GRBs are the result of delayed black hole formation, 
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which implies that the core-collapse and its subse¬ 
quent supernova may significantly precede the burst. 
The delay could be of order months to years Q, or 
perhaps as short as hours [45| . For only two of the SN 
light curves the fit indeed improved if we allowed for 
a shift in time between the onset of the burst and the 
onset of the SN (GRBs 990712, 011121). The offsets 
never exceeded 5 days, and were both negative and 
positive. However, the uncertainties in this parame¬ 
ter are large, due to the poorly sampled shape of the 
underlying supernova (e.g., |l3|). 

4. Summary and Conclusions 

Since the first clear evidence for extra light in a 
GRB afterglow light curve (GRB 980326; @|), there 
is growing evidence for several such cases. Our key 
finding is photometric evidence of a late-time bump in 
all afterglows with a redshift z <, 0.7, including those 
of the year 2003 (GRBs 030329 and 031203) and year 
2004 (GRB 041006; f42|). For larger redshifts the data 
is usually not of sufficient quality, or the SN is sim¬ 
ply too faint, in order to search for such a feature in 
the late-time afterglow light curve. This extra light 
is modeled well by a supernova component, peaking 
(l + z)(15...20) days after a burst. This, together with 
the spectral confirmation of SN light in the afterglows 
of GRB 021211, 030329, and 031203 further supports 
the view that in fact all long-duration GRBs show 
SN bumps in their late-time optical afterglows. Given 
the fact that a strong late-time bump was also found 
for XRF 030723 [To) and a less strong bump for XRF 
020903 (but with spectroscopic confirmation of under¬ 
lying SN light [dil l ) might indicate that this conclusion 
holds also for X-ray flashes (even though the finding 
of XRF-SNe might be more difficult; see jj39j). 
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